Correlation length-exponent relation for the two-dimensional random Ising model 
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We consider the two-dimensional (2d) random Ising model on a diagonal strip of the square 
lattice, where the bonds take two values, Ji > J2, with equal probability. Using an iterative 
method, based on a successive application of the star-triangle transformation, we have determined 
at the bulk critical temperature the correlation length along the strip, for different widths of the 
strip, L < 21. The ratio of the two lengths, (,l/L — A, is found to approach the universal value, 
A — 2/n for large L, independent of the dilution parameter, J1/J2. With our method we have 
demonstrated with high numerical precision, that the surface correlation function of the 2d dilute 
Ising model is self- averaging, in the critical point conformally coovariant and the corresponding 
decay exponent is n\\ =1. 
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I. INTRODUCTION 

In the presence of quenched, i.e. time independent 
disorder one generally considers different random sam- 
ples and the physical observables are characterized by 
their distribution and (n-th) moments. (For n = 1 and 
n = we have the average and typical value, respec- 
tively.) The extensive quantities, which are connected to 
the free-energy and its derivatives, have normal distribu- 
tion, thus in a single sample one measures their average 
value with probability one in the thermodynamic limit. 
These quantities are called self-averaging. There are, 
however, other observables, typically correlation func- 
tions, which are broadly distributed and the typical (or 
most probable) value is different from the average value, 
even in the thermodynamic limit. 

Such type of .phenomena takes place in disordered 
quantum systemsa, where the typical and average behav- 
ior of correlations and critical singularities are fi¥f n qual- 
itatively different- As known by exact resultspB, renor- 
malization groupu and numerical calculations!] in the in- 
finite randomness fixed point the average behaviour is 
dominated by the rare events, which occur with vanish- 
ing probability, whereas the typical behaviour is seen in 
any large sample with probability one. 

In a classical system the effect of disorder is compara- 
tively weaker, (since in quantum systems the disorder is 
strictly correlated along the (imaginary) time direction) , 
here the critical singularities are controlled by a random 
fixed point and there are usually quantitative differences 
between the average and typical behaviour. 

In this respect a well known example_is the one- 
dimensional (Id) random bond Ising modeE, defined by 
the Hamiltonian H = — JjSjSj+i. Here the spin cor- 
relation function 



j+r— 1 



G(r) = (s j+r Sj) = Yl *fc > *fc = tanh(J fc /fc B T) , 



is given as a product of random numbers and has a log- 
normal distribution. Consequently its most probable, or 
typical value G(r)t yp = exp[ln(G(r)] av = [i][ v and aver- 
age value [G(r)] av are different, even in the thermody- 
namic limit. In the following we use (. . .) and [. . .] av to 
denote thermal and disorder averaging, respectively. 

In higher dimensional classical spin systems with ran- 
dom ferromagnetic couplings, such as the random Ising 
and Q-state Potts models, the effect of disorder is ex- 
pected to be even weaker than in Id. In calculating the 
correlation function the thermal average in higher dimen- 
sions involves several random couplings, not only those 
connecting directly the two points, therefore the disor- 
der fluctuations are smoothed down. There is a class of 
random systems which, in the vicinity of their critical 
point, are homogeneous in macroscopic scales, thus the 
effect of quenched disorder is irrelevant. The correspond- 
ing criterion for weak randomness due to Harriso requires 
Q,pure w j iere a P ure i s t ne specific heat exponent of 
the pure system. 

For systems with cv purc > the disorder is a rele- 
vant perturbation so that the critical properties are con- 
trolled by a (new) random fixed point in which uncon- 
ventional scaling behavior is expected. A detailed study, 
both (field)-theoreticaB'Ll and numericafillOlIj, about 
the two-dimensional random Q > 2 state Potts model 
has revealed that the critical bulk spin correlation func- 
tion has multifractal behavior: the different moments of 
the correlation function at the critical point decay as a 
power: 



[G»] 



l/r, 



(1.2) 



with n dependent decay exponents, x^ n \ We note that 
for conventional scaling the a;^ n ^-s have no n-dependence. 

An important question concerning random mag- 
netic systems is whether the critical point correla- 
(1.2j-transform coovariantly under conformal 
ionsEJ. 



tions in Eq. ( 



(1.1) transformationsEj. Although correlations in one sample 
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are not translationally invariant the average correlations 
are translationally and rotationally invariant and - it is 
generally believed - they are also conformally coovariant. 
Indeed numerical studies in the strip and rectangle geom- 
etry for the two-dimensional Q > 2 state random Potts 
model show that average critical correlations transform 
coovariantly under conformal transformational,! Re- 
cently conformal properties of correlation functions and 
density profiles have been used tp, determine the scaling 
dimensions of different operatorsEj. 

The two-dimensional random Ising model with a purc — 
represents the marginal situation of the Harris crite- 
rion and detailed studies— have been performed to clar- 
ify its critical propertieaia. Disorder is predicted as a 
marginally irrelevant perturbation by field theory, so that 
the critical singularities of the random model are char- 
acterized by the exponents of the pure-Jsing model sup- 
plemented by logarithmic corrections^ fi^tl ., Numerical 
studies are in favour of this scenarialj'tj'EfroEj'LJ, al- 
though conflicting interpretation of the numerical results 
has also been suggestedE3o. 

Considering the spin correlation function of the ran- 
dom Ising model according to field theory the decay of 
the different moments of the bulk correlations at the crit- 
ical point are given by the power law of the pure model 
with x purc = 1/8, but the logarithmie,corrections to the 
different moments are n dependentflcJ. Numerically the 
bulk critical correlations are studied in the infinite plane 
geometry by MC simulations in Ref[ pj| , however the 
possible logarithmic corrections of the moments have not 
been analyzed. On the other hand in Ref[ [18], the 
transfer matrix method is used in the strip geometry and 
the decay exponent of the typical and average bulk cor- 
relations are deduced from the assumption of conformal 
invariance. 

In the semi-infinite geometry the critical surface corre- 
lation function has been studied in Ref [ ||l]] by the star- 
triangle (ST) method . It was found that for any dilution 
the numerically calculated average correlations are com- 
patible with the form: 

[G.WU-r-^lnr) 1 / 2 . (1.3) 

Thus the decay exponent of the critical surface magneti- 
zation is ?7|| = 1, as for the pure system, however in the 
random model there are also logarithmic corrections. 

In this paper we continue to study the surface correla- 
tion function of the 2d random Ising model. New features 
of our investigations are the following. 

i) We considered the strip geometry rather than the 
semi-infinite geometry. 

ii) As a numerical method we used an iterative proce- 
dure based on the star-triangle transformation. By this 
ST method a finite strip of random Ising model is for- 
mally transformed to a chain of Ising spins with smoothly 
inhom ogen eous bonds. Then, using the exact expression 
in Eq.(l.l), we have calculated very accurately the cor- 
relation length parallel to the strip, and studied its 
distribution and different moments. 



iii) The advantage of the ST method to the transfer 
matrix (TM) technique, applied previous for bulk cor- 
relations is two-fold. First, we could investigate larger 
widths of the strip, going up to L = 21, which is ap- 
proximately twice of the widths available by the TM 
techniques. The second advantage of the ST method 
that one can consider correlations between two largely 
separated spins, r — O(10 3 ), which is at least one order 
of magnitude larger, than for the TM method. In this 
way we obtained more accurate averaging and could go 
deeper into the asymptotic region of the correlations. 

iv) Finally, we calculated different moments of the cor- 
relation function and studied the validity of the correla- 
tion length - exponent relation, as follows from the as- 
sumption of conformal invariance. 

The structure of the paper is the following. The model 
and the ST method to calculate surface correlations are 
presented in Section 2. Our results about surface correla- 
tions and the corresponding correlation lengths are given 
in Section 3. We conclude our paper with a Discussion 
in the final Section. 



II. STAR-TRIANGLE APPROACH TO SURFACE 
CORRELATIONS 

We consider the Ising model on a diagonal strip of 
the square lattice, with i = 1,2, ...,L columns and 
j = 1,2,..., K rows. At i = 1 and i = L there 
are two (1,1) surfaces, whereas in the vertical direction 
with K L we impose periodic boundary conditions. 
The nearest neighbor spins are connected with ferromag- 
netic couplings, Ji j > 0, which could take two values, 
Ji > Ji , with equal probability. In the [thermodynamic 
limit L, K — > oo the model is self-dualcZl and the self- 
duality point: 



tanh(Ji/fc B T) = exp(-2J 2 /k B T) 



(2.1) 



corresponds to the critical point, since according to nu- 
merical studies there is one phase-transition in the sys- 
tem. The degree of dilution can be varied by changing 
the ratio of the strong and weak couplings, p = J\j J 2 . 
At p = 1 one recovers the perfect Ising model, whereas 
for p — > oo we are in the percolation limit, where T c = 0. 

For a given distribution of the couplings correlations 
between two surface spins, G s (r) = (si,j+ r si,j), can 
be conveniently calculated by the star-triangle method. 
The STjnethod was introduced by Hilhorst and van 
Leeuwenc3 and later by othersEa to study the surface crit- 
ical properties of triangular lattice Ising models with a 
layered structure. Recently the method has been gen- 
eralized for arbitrary distribution of the couplings and 
applied for the random semi-infinite Ising model in Ref[ 
pl[ , hereafter referred to as Paper I. In the following we 
recapitulate the method for the strip geometry with free 
boundary conditions at the two edges of the strip. 
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correlation function stays asymptotically invariant under 
the mapping. Then G s (r) in the original model can be 
expres sed in the form of the one-dimensional Ising model 
in Eq.(l.l) replacing Jj by the asymptotic value of the 



surface coupling Jj S \m) = Jij(m). 

For a given strip of width L the surface correla- 
tions show an asymptotic exponential decay G s (r) ~ 
exp(— t/^l), t 3> L, where the correlation length, £l, 
is approximated by 




100 



FIG. 1. Strip of triangular lattice with L = 4 layers. The 
hexagonal lattice, obtained by the ST transformation is de- 
noted by dashed lines, the triangles of the new triangular 
lattice are shown by dotted lines. In a diagonal square lattice 
the vertical couplings, which are zero in the original model, 
are generated during the iteration process. 

Central to the method is the ST transformation by 
which one replaces all right-pointing triangles of the strip 
by a star, which yields a hexagonal lattice of spins, de- 
noted by dashed lines in Fig. 1. In the second step of the 
mapping the left pointing stars of the hexagonal lattice 
are replaced by triangles resulting in a new triangular lat- 
tice, which is denoted by dotted lines in Fig. 1. Iterating 
the procedure a sequence of triangular Ising models is 
generated (m = 0, 1, 2, . . .) from the original model with 
m = 0. Neither the width of the strip nor the number of 
spins is changed under the transformation. 

As seen in Fig. 1 the surface spins of the m-th and the 
(m+ l)-th models are connected by the surface couplings 
of the intermediate hexagonal lattice and, as shown in Pa- 
per I, there is an explicit relation between the thermal av- 
erage of the surface spins in the two models. Then, from 
the fact that the surface magnetization of a finite strip 
with free boundary conditions is vanishing, follows that 
the surface couplings of the hexagonal lattice goes to zero 
as Tn — > oo. As a consequence the surface spins of the tri- 
angular lattice decouple asymptotically from the rest of 
the system. The surface correlation functions of the m-th 
and (m+l)-th triangular models are similarly related and 
one can use the results of Paper I to calculate the surface 
correlation length from this relation. One can, however, 
proceed in a simpler way noticing that the surface spin 



FIG. 2. The average surface diagonal coupling, connecting 
the surface spins to the rest of the system, as a function of 
the number of iteration, m. The calculation is performed on 
a given random sample with dilution p — 4 and on a strip of 
size L — 11 and K = 1024. Inset: The approximate inverse 
correlation length, defined in Eq.(2.2), as a function of the 
iteration parameter. 
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(2.2) 



and limm._i.oo ^l(to) = £l. Averaging over different disor- 
der realizations one obtains the typical correlation length: 



c*yp 



L I av 



(2.3) 



In physical applications one should average the n-th 
power of the correlation function the asymptotic behavior 
of which, [G s {r) n ]\i n ~ exp(— r defines the corre- 
sponding correlation length, £^ . As already mentioned 
for n — 1 and n = we obtain the average and typical 
correlation lengths, respectively. From the different mo- 
ments of the distribution of the inverse correlation length 
p(1/£l) on e obtains £r in a cumulant expansion: 
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(2.4) 



For isotropic systems the correlation length parallel to 
the strip, £^ and the width of the system, L, are asymp- 
totically proportional and for conformally-invariant sys- 
tems their ratio takes the universal valuet_: 
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(2.5) 



Here Xg is the anomalous dimension of the surface mag- 
netization, defined through the asymptotic decay of the 
critical surface correlation in the semi-infinite geometry, 

Thus Xg is the surface counter- 



,1/n 



-2x\ 



[G s (r)"] av 

part of x^ in Eq.(1.2) and satisfies the scaling relation 
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III. RESULTS 

We studied the spin correlations of the random Ising 
model on the (1,1) surface of the square lattice by the 
ST method. Evidently the original model with to = 
can be considered as a special triangular lattice model 
with vanishing vertical bonds. During iteration, however, 
non-zero vertical couplings are generated so that also the 
surface couplings (to) become non-zero. In the ac- 
tual calculations we considered strips of width L = 21 + 1 
up to L = 2lEJ, whereas for the length of the strip, K, 
the condition K ^> L is always satisfied. Typically we 
took K = 1024 and checked that the numerical results 
are insensitive on the variation of K in this region. 
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FIG. 3. Probability distribution of the inverse of the corre- 
lation length for L = 11 and for different lengths of the strip, 
K. We studied 20000 samples for each lengths at a dilution 
p = 4. The full lines represent the Gaussian approximations 
to the distributions. 

As indicated in the previous chapter, under iteration 
the surface spins asymptotically decouple from the rest 
of the system and the surface correlations have one- 
dimensional character. For an illustration we have cal- 
culated the average value of the first diagonal coupling, 
[j( d )(m)] av , connecting the first and second line of spins, 



as a function of the iteration, m. It is given in Fig. 2 
together with the correlation length, £,L( m ), as defined 
in Eq. ( |2.2| ) between two spins of maximal distance, 
r = K/2. As seen in Fig. 2 both [J (d) (m)] av and £l(to) 
approach their limiting values rapidly, exponentially with 
to. Analyzing the iteration equations in Paper I one can 
show that to ~ L 2 iteration steps are needed to reach the 
asymptotic region, which is indeed verified numerically. 

For a given dilution, p, the correlation length, 
shows variation from sample to sample. The distribu- 
tion of the inverse correlation length, p(1/£l), obtained 
over 20000 samples is shown in Fig. 3 for different lengths 
of the strip, K . As seen in Fig. 3 the average v alue of 
£l, defining the typical correlation length in Eq.(2.3) is 
independent of K, whereas the width of the distribution 
is decreasing with the length of the strip as 1/ y~K. This 
observation is in agreement with the cumulant expansion 



in Eq.(2.4) and with the fact that £^ is asymptotically 



independent of K. The distribution p(l/£i) is found ap- 
proximately Gaussian, however for finite strips there is 
always some deviation from the normal distribution. 
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FIG. 4. Ratio of the typical correlation length, £/ p , and 
the with of the strip, measured in units I = (L — l)/2 for 
different dilutions. The dashed lines are guide to the eyes, 
representing the finite-size corrections as 0(1/1). 

Next, we are going to study the L dependence of the 
typical correlation length, £^ yp , for different values of the 
dilution p = 1, 2, 4 and 10. First we note that in order to 
use the same lattice units in the vertical and horizontal 
directions one should replace L = 21 + 1 by I. In Fig. 4 
the ratio (tt/2)^ p /I is plotted against l/l. As seen in 
Fig. 4 in the range of dilution we worked, the correlation 
length is monotonically decreasing with p, whereas there 
is an approximate linear 1 jl correction to the ratio for all 
values of p. The asymptotic value of the ratio is found 
dilution independent, we estimated as 
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(3.1) 



IV. DISCUSSION 



We have also studied the different moments of the cor- 
relation function and calculated the corresponding av- 
erage correlation length, £^ . The ratio (n /2)^£' /I is 
found to have a strong size dependence, much stronger 
than for the typical correlation length. In this case the 
finite-size corrections are approximately logarithmic, as 
can be seen in Fig. 5. For any finite L the corrections are 
increasing with the moment, n, however the asymptotic 
value of the ratio tends to the same universal value 



lim — 

L — >oo 2 
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= 1±0.03 



(3.2) 



as for the typical ratio in Eq.( 3.1). We note that the 
error of the estimate in Eq^jB^) is increasing with n. 
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FIG. 5. Ratio ol the correlation length corresponding to 
the n-th moment of the correlation function, £^ , and the 
width of the strip in units of I = (L — 1)/2 at a dilution p = 4. 
The dashed line is guide to the eyes, representing logarithmic 
finite-size corrections as 0(1/ In/). 



At this point we compare our results in Eqs. (3.1) and 
(|3.2[) to that obtained in the semi-infinite geometry in 



Eq.(1.3) and comment about the val idit y of the correla- 



tion length-exponent rel ation i n E q.(2.5). The universal 
ratio found in Eqs. (3d) and (3^) according to the cor- 
relation length-exponent relation corresponds to a typ- 
ical and average surface magnetization decay exponent 
of ?7|| = 1 in complete agree ment with the result in the 
semi- infinite geometry in Eq.(1.3). Furthermore, for the 
average correlations we observed logarithmic corrections 
in both geometries. Thus we can conclude that the cor- 
relation length-exponent relation is valid for the typical 
and average surface correlations of the random 2d Ising 
model, thus they transform coovariantly under conformal 
transformations. 



In the previous section the surface correlation function 
of the 2d random Ising model is studied in the strip ge- 
ometry. In particular we have calculated the ratio of the 
correlation length, £ L , obtained from the average of the 
n-th moment of the surface correlation function, and the 
width, L, of the system. We found that asymptotically 
this ratio goes to a universal value, irrespective of the de- 
gree of dilution, p, and the value of the moment, n. For 
typical correlations the correction terms are of 0(1/L), 
whereas for the average, n = 1, and for the higher mo- 
ments, n > 1. The finite-size corrections are logarithmic, 
they are in the form of 0(1/ In L). 

Several qualitative features of the above results can 
be understood by analyzing the ST iteration procedure. 
As mentioned before the surface spins are asymptotically 
decoupled after m ~ L 2 iteration steps, when the expres- 
sion of the new surface couplings, J^(m), contains a 
set of the original couplings, J,-j, taken from a region 
of 1 < i < L and k — L < j < k + L. Consequently 
the jjf\m) are correlated for short distances, but they 
are practically independent between two sites which are 
separated by at least a distance of 0(L). Now we can 
use an approximate coarse-grained description: introduce 
block-spin variables to replace each 0(L) number of sur- 
face spins which are connected by correlated couplings. 
The new couplings between the block-spins, J k , , are ap- 
proximately independent random variables and they sat- 

isfy tanh( J^, /fcgT) = 0(1), since with the original cou- 
plings tanh( Jl s) /k B T) = 0(1/L). Then from Eq.(||) 
follows that the correlation length of the system isself- 
averaging, in any sample the expression in Eq.( [2.2| ) in 
the thermodynamic limit goes to the average correlation 
length, [^ijav, with probability one. Evaluating the vari- 
ance of 1/£l in the coarse-grained picture one gets: 



typ 



rL 



(4.1) 



thus the first correction term in the cumulant expansion 
in Eq.(2.4) is of 0(1/L), as it should be. The eventual 
correlations between the block-s pin couplings will tend 
to reduce the value of j3 in Eq.(4d). According to our 
numerical studies ~ 1/ In L and this effect is the source 
of the logarithmic corrections in the random Ising model. 

The results in the semi-infinite geometry and in the 
strip geometry are in complete c orre spondence. Compar- 
ing the conformal result in Eq.( |2.5| ) with the numerical 
estimates in Eqs. (3.1) and ( |3~2] ) one obtains the following 
conclusions. 

i) The correlation length-exponent relation is valid for 
the random Ising model, thus the (surface) correlations 
of the system are conformally coovariant. 
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ii) For the typical and average correlations at the criti- 
cal point the decay is given by the same exponent, which 
does not depend on the degree of dilution. Consequently 
there is no multifractal behavior for the critical correla- 
tions of the model. 

iii) The typical surface correlations are free of logarith- 
mic corrections, whereas the average correlations and the 
higher moments are subject of logarithmic corrections, 
the strength of those is increased with n. 

iv) Finally, our numerical results give strong and accu- 
rate numerical support to the field-theoretical conjecture 
that the random and pure Ising models in 2d belong to 
the same universality class. 
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